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We predict a novel effect in a quantum two-level system (TLS) coupled to a resonant cavity. By 
bringing the TLS in and out of resonance with the cavity by a series of N rectangular bias pulses 
(the length of the mth pulse scaling as 1/^/m), we will coherently excite the JV-photon state, TV ) , 
of the cavity only if the TLS was initially in an appropriate quantum state ("go" state). Otherwise 
the number of photons in the cavity will remain small compared to N {selective amplification). If 
the TLS was in a coherent superposition of the "go" and "no go" states, the cavity will be in a 
superposition of states, in which the state |JV } will enter with the same weight as the initial "go" 
component. The effect is due to \/N + 1-dependence of the Rabi oscillation frequency on the number 
N of photons in the cavity. It is stable with respect to noise, pulse shape, finite temperature, TLS 
decoherence and TLS detuning from resonance with the cavity. The effect can be used as a means 
to read out a quantum state of a qubit coupled to a resonator. 

PACS numbers: 



INTRODUCTION 

A two-level system (TLS) coupled to an oscillator is a 
textbook example in quantum optics [1], where dressed 
states, Rabi oscillations and entanglement between an 
atomic degree of freedom and cavity modes can be in- 
vestigated [2]. More recently, this model was used to 
describe the behaviour of superconducting quantum bits 
interacting through some kind of a resonator or a bus 
[3-7] . Drawing the analogy with cavity QED [2] , one can 
entangle qubits with each other by consecutive coupling 
them to the resonator, creating in the process a coherent 
superposition of number states in the resonator, e.g. \n ) 
and \m ), with close values of n, m. 

If we could put the resonator in a superposition state 
a\N ) + (3\M ) , with a and (3 depending on the state of 
a TLS coupled to it, the latter can be determined by 
measuring the resonator, if only n is above the detection 
threshold, m is below it, and the difference \m — n\ ex- 
ceeds the noise level. We show that such a phenomenon 
actually takes place, under the condition, that the en- 
ergy levels of the TLS can be made (almost) degenerate, 
with tunneling splitting A <C luq, the frequency of the 
resonator [*1. 



The Hamiltonian of the system is 

H = uj G a)a + H q + H u 



Here the TLS Hamiltonian is 



H a 



1 



1 



:£<7, 



Aa x . 



(1) 



(2) 



It is written in the basis of "physical" states, which are 
the eigenstates of the qubit at infinite bias e. (For a 
superconducting flux qubit [10], which is a loop with 
a circulating supercurrent, such states correspond to 
clockwise and counterclockwise current directions. For 
a charge qubit [11], those are states of a small supercon- 
ducting grain with a definite number of Cooper pairs.) 
We will denote the corresponding states of the TLS- 
resonator system by \L, n ) and \R, n ) . In the same basis, 
the interaction is 



Hmt = (A z <r 2 + \ x o- x )(a f + a). 



(3) 



In the absence of small parameters, one would be re- 
duced to solving numerically the equation of motion for 
the density matrix of the system, 



ip = [H, p] 



(4) 



This amplification scheme for detection of a qubit state was first 
suggested by the authors in Ref. 8. Similar technique was pro- 



posed in Ref. 9 as a means of creating arbitrary Fock states in a 
microcavity. 




FIG. 1: Physical picture of selective amplification 
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In the near-resonant regime, f2(e) = VA 2 + e 2 w 
wo, the system is described by the Jaynes-Cummings 
Hamiltonian[l] (Fig.la,d,e): 



■yjRWA 
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(6) 



(7) 

Here the Pauli matrices a x , a z , a± — a x =F icr^ are in the 
eigenbasis of biased qubit (so that the states of the system 
are now |0, n ) and |1, n ) ). The subscripts remind, that 
the basis is bias-dependent; obviously, we can make a 
convention that |0,n } = \R,n ) , |l,n ) = |£, n ) (and 
the opposite when e = -co). The effective coupling is 



9(e) 



A Z A + X x e 
fi(e) ■ 



(8) 



The block-diagonal matrix V(s) = ©v(e) realizes the 
transformation between this "energy" and the initial 
"physical" basis: 



v(e) = 



= J_ / y /1 + £/»(£) - Vl-g/n(e ) 



^2\^/l-e/n(e) y/l + e/n(e) 



(9) 



(for the time intervals less than the dephasing and decay 
time in the system, which we assume to be much longer 
than other time scales, Ti/uiQ, 7i/A). We will later present 
the numerical solution of this equation in the presence of 
noise, to demonstrate the stability of the effect. But first 
let us use the rotating wave approximation (RWA) [1], 
which provides a clear physical picture of the process and 
turns out to be quantitatively accurate for the relevant 
choice of parameters. 



The Hamiltonian i-[^ WA can be diagonalized (e.g. [1]). 
Its eigenstates are superpositions of states with the same 
excitation number (e.g., |0, n ) e and |l,n — 1 ) E ), and if 
initialized in the state with a fixed photon number in the 
resonator, the system will undergo quantum beats with 
frequency, dependent on n, which can be considered as 
Rabi oscillations of the TLS in the presence of the res- 
onator field (see Fig. la,d,e). To be precise, the evolution 
operator, exp[— itH^ WA ], up to a phase factor, is given 
by a block-diagonal matrix 



SELECTIVE AMPLIFICATION IN THE 
ROTATING WAVE APPROXIMATION 

The condition A <C wo allows us to consider two sepa- 
rate RWA regimes. In the small bias regime, e <C uo, the 
TLS is effectively decoupled from the resonator due to the 
large discrepancy in the relevant frequencies (Fig.lc): 



/fflq>[-itf(e)] 



S(e,t) = 



So(e,t) 



si(e,t) 



V 
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(10) 



Here 6(e) = lo$ — 0(e) is the detuning between the TLS 
and the resonator. The 2x2 matrices Sj,j — 0, 1, . . .: 



J 



a j = e- ijUo *t 
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( QQ,e)t 
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iSM „: n ( QU,e)t \ -2ig(e)VJfT . ( Q(J,s)t 
Q(j,e) Sm \ 2 J Q(j,e) &H1 \ 2 



S|Q y^^-^j 



( QU,e)t " 



cos(^i 



Q(j» Sm ^ 2 



i = ,ju t( Cj (*) - ify (*) Vin 
'-' ! -^(t)c*(t) / iiJ 
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describe the quantum beats between the states \l,j 



and |0,j + 1 ) , belonging to the same (j + l)th sector 
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FIG. 2: Probability distribution of occupation numbers of the 
resonator, P n (m) = p T nn (m), as a function of number of bias 
cycles, m, calculated using RWA. In the units of A, the pa- 
rameters are lu — 10 (resonator frequency), A z = 0, X x — 0.8 
(TLS-resonator couplings), 5 = (perfect resonance), |e| = 
9.9499. The length of the mth bias pulse is 3To/y/m, where 
To = 5.9202 is the Rabi half-period for empty resonator, (a) 
Zero temperature: empty resonator; initial density matrix of 
the TLS is p TLS (0) = |1 ) e ( l| e ("go" state, left), |0 ) e ( 0| e 
("no go" state,right). (b) Same for finite temperature (res- 
onator in equilibrium at T — 1); (c) Same for T = 20. 



(labeled by the total number of excitations in the sys- 
tem). The corresponding frequency is Q(j,e)/2: 



Q(j,e) = y/P(e)+4g*(e)(J + l). 



(12) 



The vacuum state, |0, ) e , is in a sector of its own and 
acquires only a phase factor. 

Now we can discuss qualitatively the evolution of TLS 
and resonator. Suppose the system is initially in a pure 

state 

oo 1 
i— a— a,j 

and assume for simplicity exact resonance, 5(e) = 0. It 
is convenient to introduce complex two-dimensional vec- 
tors Cj = (Cqj, Cij-i). Then, according to (11), after 
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FIG. 3: Importance of phase randomization: time dependence 
of P„ for the same parameters as in Fig. 2(a), but with the bias 
pulse length T m — To/^/m. The "no go" state is also amplified 
and can not be reliably distinguished from the "go" state. 
Further increase of the /^-factor beyond 3 will not qualitatively 
change the outcome (not shown). 
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FIG. 4: Role of detuning: w = 10, \ z =G,\ X = 0.8, 5 = 1. 



time T p — tt/Q(p — l,e) the vector C p flips (up to a to- 
tal phase factor), while all other vectors C^ p turn by 
some angle determined by the ratio T p /T[. If wait for a 
large odd multiple of T p , the rotation angles of all / ^ p- 
vectors become quasirandom. This will be important for 
the following. Still, the maximum number of photons in 
the reservoir cannot change by more than one. (In the 
simplest case, we transformed |l,p— 1 ) — ► \0,p ).) 

In other words, we coherently discharged the TLS into 
the resonator (Fig. lb). In order to pump more pho- 
tons, we must recharge the TLS by switching between 
the excitation-number sectors (e.g. \0,p) — > \l,p )). 

There are two possibilities. We can suddenly lift the 
bias, £ — > 0, decoupling the resonator from the TLS 
(Fig. lc) (as suggested in Refs. 8, 9). The TLS state 
|0)~|i?)isa superposition of the eigenstates of unbi- 
ased TLS, |0(1) ) = QR ) + (-)|£ ))/\/2. After an odd 
multiple of To = 7r/A, it will approximately become the 
\L )-state, and if then we reapply the same bias e, the sys- 
tem will find itself in the desired state | l,p ) (Fig. Id). 
The accuracy of this transformation is easily calculated. 

Another method is instead to reverse the sign of the 




bias [f], Fig. Id. Now we simply switch between the 
e-basis and (— e)-basis, with the transformation matrix 
V(-e)Vt(e) = U{e) = u(e), where 



(13) 



if A <C e, 0(e), and come to the "recharged" state 

|l,P)-e- 

We will be discussing the later method, as more 
straightforward. Nevertheless, if the actual physical sys- 
tem described by a TLS does not allow to switch the sign 
of the bias, but only to bring the levels close to the de- 
generacy, the former method can be applied to the same 
effect. 

Note that in both cases the TLS-switching transforma- 
tion does not depend on the resonator occupation num- 
ber. On the contrary, the transition process of the addi- 
tional photon to the resonator is sensitive to it; therefore 
we need now to wait an odd multiple of T p+ \ to flip the 
vector Cp+i; all other vectors will again turn by quasir- 
andom angles. 

Now we see how the selective amplification 
works. If apply to the TLS a series of res- 

onant bias pulses (e, — e,e, — e, . . .) of duration 
[HiT p , n^Tp+i, ^T p+ 2, ^T p+ 4, ...) (fij is an arbi- 
trary odd number, then the weight of the |l,p)- 
componcnt of the initial wave function of the system 
will be transferred upwards (up to a phase factor): 
C\p — > Co, p +i — > Ci tP +i — > Co, p +2 — ► • • ■ On the 
other hand, all other weights will be quasirandomly 
redistributed leading to, at best, diffusive growth of 
the maximal occupation number of the resonator. If 
we measure this number by a detector with ./V-photon 
threshold after N — p + s steps, the probability of nonzero 
reading will be given by the initial weight of the state 
\l,p ) , as long as the dispersion of the average photon 
number is less than s 2 . Therefore we can measure the 
initial quantum state of the system. 

The effects of finite temperature can be easily taken 
into account under a natural assumption that the energy 
relaxation time for our system is much longer than the 
observation time. (Otherwise it would be impossible to 
observe quantum coherence at all.) The system will be 
then described by the density matrix pit), with the initial 
value p(0) = /9 TLS (0) ® p r T , where pi T is the equilib- 
rium density matrix of the resonator at temperature T. 
Its time evolution is given by 



(ip(^-^-i))p(o)(n 5 ( 



l &i i vi £?' 



4 



(14) 



The results are presented in Figs. (1-6), where the di- 
agonal elements of the resonator density matrix after 
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FIG. 5: Selective amplification of a superposition state, 
p TLS (0) = |(|1 ) + |0 » e (( 1| + < 0|) e . The parameters are 
ui = 5, X z = O.A^ = 0.1, 5 = 0, \e\ = 4.8990. The length of 
the mth bias pulse is 3To/\/m, where To = 48.0956 is the Rabi 
half-period for empty resonator. Temperature T — 0; T — 1 
(left), T = 10 (right) (the pictures for T = and T = 1 are 
visually indistinguishable on the chosen scale -cf. Fig.2(a,b)). 
Note that the effect survives in spite of smaller coupling and 
smaller mismatch between u) and A. The state of the res- 
onator after 20 pulses (rightmost column) is approximately 
(|0) + |20))/v / 2. 
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[f] AZ is grateful to B. Wilson for pointing it out. 



FIG. 6: Role of decoherence (the state of TLS is measured 
after every bias cycle, as described in the text): u = 10, 
A z = 0,Ax = 0.8, 5 = 0. The temperature T = (top), 
T = 20 (bottom). 



the mth bias step, P n (m) = p r nn (m), are plotted. Here 
p r (t) = tvTLsp(t)- A runaway increase of the resonator 
occupation number for the "go" state is clearly seen. 
(The detailed description is given in the captions.) 

The calculations were performed for different choices 
of parameters. In particular, the choice for Fig. 5 is com- 
parable to the parameters of a charge qubit coupled to a 
strip line [7]. 

Fig. 3 demonstrates the crucial importance of phase 



randomization: after only half Rabi period, the "no go" 
state is not sufficiently wiped out by destructive interfer- 
ence and also produces a significant runaway component. 

On the other hand, the stability of the effect with re- 
spect to finite temperature, detuning, and variation of 
system's parameters is remarkable. Note that the super- 
position initial state of the TLS, with the vacuum initial 
state of the resonator (Fig. 5), leads after N bias pulses 
to the resonator state of the form a\N ) + (3\M ) plus 
small admixture of other low-lying states. 

The role of decoherence in the TLS warrants special 
attention. Each bias application increases the risk of de- 
stroying its phase coherence. In order to model this pro- 
cess, we repeated the above calculations, wiping out the 
elements of p corresponding to coherent TLS states (in 
our convention, all pij with i and j of different parity) af- 
ter each bias application. The rest off-diagonal elements 
of p correspond to quantum correlations in the resonator, 
which will likely have a much longer decoherence time 
than the TLS [5, 12]. The results presented in Fig. 6 
show that there is virtually no effect of such decoherence 
on the evolution of the "go" state, while the structure of 
the "no go" probability amplitude is smeared out. This 
is physically clear: after a bias cycle, the "go" state is 
factorized, and measuring the TLS state (which was es- 
sentially done to model decoherence) does not affect the 
resonator state, which keeps all the relevant quantum cor- 
relations. For the "wrong" states the process reduces to 
an additional averaging, wiping out the remaining struc- 
ture of the density matrix. 



ANALYTICAL SOLUTION 

The decoherence model used in the previous section 
provides a crucial simplification, which allows an analyt- 
ical solution of the problem. 

Let us introduce the rotated density matrix of the sys- 
tem after the nth step: 

p{t n ) = V{e n+1 )p{t n )V\e n+1 ). (15) 

Substituting here the time evolution (14) of p, we obtain 

p(t n ) = V(£„+i)V t (e„)5(e„, t n - £„_i) x (16) 

xp(t n )S f (e n ,t n - i„_i)V(£„)V t (e n+ i). 

In order to proceed, we must make some simplifications. 
First, assume that e„ = — e n —i and |e| 3> A. This 
replaces the products V(t n )V^(t n -i) with simple block- 
diagonal icr z -matrices. Second, use RWA. Third, assume 
that the TLS loses coherence after every cycle, i.e. keep 
only the elements of the desnity matrix with the indices 
of the same parity. Then we obtain for the relevant ele- 
ments of the matrix p (/z, A = 1, 2, . . .: 



[/5(£n+l)] 2M ,2A - 



+ [/5(*n)]2„+l,2A+l Vl(*«)SA-l(*n)e- i( ' i - AV0(t ''- t "- l) ; 

[P(*n+1)]2 M +1,2A+1 = 
= [p(in)] 2/i+ l,2A + lC;-l(*n)C A _ 1 (t„)e- i (' i - A )^^- t «- 1 ) + 

while [ J 5(t n +i)] 11 = [/5(£ n +i)] n . Here Cj,Sj are the ele- 
ments of the evolution matrix (11). Different diagonals of 
the density matrix thus transform independently We are 
interested primarily in the main diagonal, the elements of 
which, 7VJ = Pjj(t n ), describe the photon numbers in the 
resonator (the probability for the resonator to contain p 
photons at the moment t n is P n {p) = 7?-2p+i + 7^2p+2, 
p = 0, 1, . . .). The recursive relations for IZ's are 



f K 



71+1 



\c (t n )\ 2 n% + s (t n ) 2 n%; 
ft?; 

(p=l,2,...) 



ft™ +1 



K. 



n+l 
2p+l 



(17) 



|Cp(*ra)| 7c 2p _|_2 + O p {t n ) TZ-2p+3i 



, ft 2 p+2 — l^p-l^n)! 7?-2p+l + dp-l{t n ) ft^p! 



They can be further simplified if take into ac- 
count the quasirandomization. Indeed, the coefficients 
\C p (t n )\ 2 , S p (t n ) 2 contain cos 2 (sin 2 ) (Q p (t n )t n /2) , and 
unless n = £>[$], the oscillating terms will average to zero: 



<|Cp(U 2 ) = ^(i + gJ-^)(i-<V0 + 



(18) 



1 s 2 

\(H 2 P )(i- 



( s P (t n f ) 



14{p+l)g(t r 



!p\tn) 



Jpm 



1 4(p+l)g(t n ) 2 

4 (°cp\J'n) 



( 1 - Spn) + (19) 



1-i 

4 



~>pm 



where we introduced the noise through the dispersion of 
phase gained by the "go" state on the nth step,( 5<p 2 ). 
The above expressions are explicitly unitary. 

Now for illustration assume (5 = (exact tuning) and 
( S<\> 2 ) = (no external noise). Now the coefficients 
enp = ( \C P (t n )\ 2 > = (l/2)(l-S pn ), s np = (( S p {t n ) 2 > = 
(1/2)(1 + Spn), and defining c„,_i = 1, s n -\ = 0, we 
simplify the recursion (17) to 



-TpTl + l 

/ ' c 2p+l 



Cnp '^-2p+2 



S r jp/C 2 p_|_3; 



„ , T? n I „ -p n . 

2p+2 — c ",p-l' v 2p+l ~ s n,p-l' v 2pi 

(n,p = 0,1,...). 



(20) 



First, note that there is a solution T^p+i 



0,K. 



2p+2 



It describes the amplification of the " go" state: the 



[J] Or, more generally, n = p — q. This would correspond to a pulse 
sequence, which amplifies the initial state with q photons in the 
resonator. 



distribution of photon number in the resonator P n (p) = 
5 p ,n- We neglected noise and detuning, therefore we did 
not reproduce the slow probability spread seen in numer- 
ics (cf. Fig. 2); the average photon number after n steps 

Second, for "no go" terms, n ^ p,p—l, Eq.(20) reduces 



to 



qpn+1 ll?'- 

/-c 2p+l — 2 /v 2 



2 ,v -2p+2 
1-pn 
2p+2 — 2'^2p+l "T" 2 /v 2p' 



•2p+ 



2 '^2p+3' 
I V n 



(21) 



which yields in the continuous limit the diffusion equation 
for the photon number distribution function: 




d n p n ( P ) = -#Lp„(p), 



(22) 



with the reflective boundary condition at p = 0. The 
solution for the initial condition S(p — po), 



Pn{p\Po) 



Op— po) 



/-v/in, 



(23) 



describes the diffusive spread of " no go" states seen in the 
numerics. In particular, for po = the average photon 
number after n steps ( p ) n n ° s ° = s/n/ir. 



EXACT NUMERICAL RESULTS AND NOISE 
EFFECTS 

A direct numerical solution of the equation (4) is also 
possible. In Fig. 7 we show the resulting time evolu- 
tion of the average occupation number of the resonator, 
( N ), for "go" (curve 1) and "no go" (curve 2) states at 
zero temperature. Unlike RWA, here we could use the 
realistic, smooth pulse shape. (Fig. 8). Still, the results 
obtained in RWA (squares) are remarkably accurate. 

So far we did not discuss the effects of external noise. 
The most obviously dangerous will be the fluctuations of 
amplitude and duration of the bias pulses. Their action 
is similar, since either primarily introduces the phase er- 
ror between the Rabi oscillations of the TLS and applied 
pulse sequence. We therefore investigate the case of ran- 
dom pulse duration. Different models can be considered 
depending on the expected source of fluctuations (e.g. 
telegraph noise, 1/f noise etc). For our purposes it is 
enough to consider here a Gaussian noise with standard 
deviation St, most likely determined by the properties of 
the pulse generator. 

The influence of such noise is presented in in Fig. 9, 
where both ( N ) and the r.m.s. deviation SN are plotted 
as functions of time for different values of St. The effect 
is again stable. An obvious estimate for the maximal 
number of excitations, which can be coherently pumped 
into the resonator, is 



N„ 



(n/sty 



(24) 



FIG. 7: Selective amplification: exact numerical results for 
u = 5,X x = 0,X x = 0.1,5 = 0,T = 0, for "go" (curve 1) and 
"no go" (curve 2) states. The average resonator occupation 
number ( N ) is plotted versus dimensionless time r = h/A. 
Squares show the results obtained in RWA. 
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FIG. 8: The bias pulse sequence used in the calculations of 
Fig.7. The (absolute) sharpness of individual pulses is inde- 
pendent on the pulse duration. 




FIG. 9: Effect of noise on selective amplification. Equation 
(4) was solved with the pulse sequence of Fig. 8, modified by 
the Gaussian fluctuations of the pulse duration with standard 
deviation St (in units h/A). The average occupation number, 
( N }, and r.m.s. deviation, SN, are plotted for St — (curves 
1 and 2), St — 0.01 (curves 3 and 4), and St — 0.1 (curves 5 
and 6 respectively). 



After that the phase error introduced by the noise will 
be of the same order as the phase gain during the cycle, 
thus destroying the constructive interference required for 
the effect. 

In conclusion, we predict a new effect, selective ampli- 
fication of a quantum state of a TLS, allowing a coherent 
pumping of excitations in a resonator by a special se- 
quence of bias pulses applied to the TLS. The effect can 
be used in order to read out a state of a quantum system 
coupled to a resonator. The effect will be stable with 
respect to external noise, finite temperature, detuning 
between the TLS and resonator frequency, and decohcr- 
ence in TLS. 
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